Angular distributions in three-body baryonic B decays 
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We study the angular distributions of the baryon-antibaryon low-mass enhancements in the three- 
body baryonic B decays of B — > ppM (M = K and tt) in the framework of the perturbative QCD. 
By writing the most general forms for the transition form factors of B — > pp, we find that the angular 
distribution asymmetry in ppK~ measured by the Belle Collaboration can be explained. We 

give a quantitative description on the Dalitz plot asymmetry in _B~ — > ppK~ shown by the BaBar 
Collaboration and demonstrate that it is equivalent to the angular asymmetry. In addition, we 
present our results on ppKs and B~ ppTv~ and we obtain that their angular asymmetries 

are —0.35 ± 0.11 and 0.45 ± 0.10, respectively, which can be tested by the ongoing experiments at 
the B factories. 
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The three-body baryonic B decays oi B ^ ppM [M = 
K'^*\tt) with a near threshold enhancement in the pp 
invariant mass spectrum have been observed by the 
Belle [l|, and BaBar fsl Collaborations. In fact, the 
dibaryon threshold enhancement has been also observed 
in other decays, such as the charmless modes B^ 
Ap7r+, B- Apj and B~ AAJsT" 0, ^ and the 
charmful ones B° — > ppD^*^", — > npD*^ and B~ 
A+p7r~ [1] as well as the J/^ decays J/^E* pjr/ and 
J/* ^ ApK- 0]. Theoretically, the dibaryon threshold 
enhancements in the three-body baryonic B decays were 
first conjectured in Ref. Q. Subsequently, various inter- 
pretations, including models with a baryon-antibaryon 
bound state or baryonium exotic glueball states 
IC, 11 1, fragmentation [llj and final state interactions 
12l | have been proposed. Furthermore, the enhancements 
in B ^ BB'M (B,B' = p, A and M = DJ<,Trrf) have 
been successfully understood l^, Q M, \m m within 
the framework of the perturbative QCD (PQCD) based 
on the QCD counting rules l^, 3] due to the power ex- 
panded baryonic form factors. On the other hand, the 
threshold enhancements in the J/^* decays are still in 
favor of some bound states . 

To find out the origin of the threshold enhancements 
and distinguish among the above theoretical models, var- 
ious experimental studies have been performed (20j . The 
Belle Collaboration has studied the angular distribution 
asymmetry Q in the helicity frame for the decays of 
B~ PpK~, B° ppKs and 5° Ap7r+, while the 
BaBar Collaboration has measured the Dalitz plot asym- 
metry Q in the decay of B^ — > ppK~ , which are sensi- 
tive to the physical nature at the threshold as well as the 
decay mechanism. In particular, both data support the 
quark fragmentation mechanism but disfavor the gluonic 
picture. So far, there is no consistent understanding of 
the two asymmetries and thus one cannot directly com- 
pare the data between the two Collaborations. 

The angular distribution in i?" —f Apir^ has been ex- 
amined in Refs. Qlill based on the QCD counting rules 



and the result is consistent with the data [3]. However, 
such an approach faces a difficulty [21], |22] to understand 
the shapes for the modes of i? — > ppK due to the 
lack of understanding of the baryonic transition form fac- 
tors in _B pp. Moreover, quantitative descriptions of 
the Dalitz plot asymmetries are not yet available besides 
qualitative quark fragmentation mechanism fll|. 

In this letter, we will concentrate on the three-body 
charmless baryonic decays oi B ^ ppM {M ~ K,tt). 
Wc will first write down the most general form of the 
B —>■ BB' transition matrix element and then fix the 
unknown form factors by the experimental data. In our 
study, we will use the QCD counting rules as well as 
the SU{3) flavor symmetry to get and relate the behav- 
iors of the form factors. Wc will demonstrate that when 
these form factors are constructed, the data of the an- 
gular distribution and Dalitz plot ^ asymmetries in 
B~ — > ppK^ can be understood. We will show that the 
two asymmetries are equivalent by describing the same 
physics. Furthermore, we will extend our investigation 
to the decays of B^ PpKs and B~ ppn~ . 

From the effective Hamiltonian at quark level [2^ , the 
decay amplitude of B~ ppK~ is separated into two 
parts, given by 
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where G_f is the Fermi constant, Vq-q. are the CKM ma- 
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terms of the effective Wilson coefficients c^^'^ , defined in 
Refs. [H, [iil, and the color number Nc- In Eq. H]), we 
have assumed the factorization approximation. As seen 
in Eq. ([T]), the current part of the amplitude involves 
time-like baryonic form factors from {pp\{qq)v,A\0) a-nd 
B —f K transition form factors, which can be referred 
in Refs. 17| and (2^, respectively. On the other hand, 
for the transition part of the amplitude, we need to eval- 
uate the transition matrix elements of i? — > pp. Based 
on Lorentz invariancc. the most general forms for the 
B —!■ BB' transition matrix elements due to scalar, pseu- 
doscalar. vector and axial- vector currents are given by 



(BE'l^'-li?) = *S(pb)[/a + /p]75«(Pb') : 
{BB'\P>'\B) = tuipTi)[fv ^+fsHPB') , 
{BB'\V^\B) = zii(pB)[5i7M + 92i(Jt,vP'' + gzPf, 

+54(Pb' +Pb)m +55bB' -Pb)/.]75w(Pb') 1 
{BB'\Al\B) = zm(pb)[./i7m + /2*fTM-p" + hp^. 

+fi{Pii' +Pb)m + /5(Pb' -Pb)m]w(Pb') : 



(2) 



respectively, where S'^ = qb^ P'^ = qjsb, ~ (77^6 and 
= 97^756 with q = u,d and s, p ~ pb^Pb —Pb' is the 
emitted four-momentum, and gi and fi are the the form 
factors to be determined. Here the parity conservation 
in strong interactions is used. We note that the forms 
for the scalar and pseudoscalar currents in Eq. ([2]) have 
been studied (e.g., see Ref. [S]), but those involving the 
vector and axial-vector currents in Eq. ([2]) have not been 
given in the literature previously. 

At the scale of to;, ~ 4 GcV, the t = {pb +Pb' )'^ depen- 
dences of the form factors in Eq. ^ can be parametrized 
according to the power counting rules of the PQCD 
26l . 27 1 based on the hard gluons needed in the process. 
For example, in the transition of i? ^ BB', three hard 
gluons are needed to produce BB', in which two of them 
create the valence quark pairs and the third one is re- 
sponsible for kicking the spectator quark in B [2l| . Since 



each of the gluons has a propagator ~ 1 /q with q pro- 
portional to the momentum of the BB' pair, all of the 
form factors /; and gi have to fall off as 1 /t^ , in which the 
propagator can be realized to contain a zero-mass pole 
inducing threshold enhancement Q . Explicitly, we have 



gi 



i3 ' ' t3 



(3) 



where Cg^ and Cf- are new sets of form factors to be 
determined. 

From equation of motion, we can relate the form fac- 
tors in Eq. ([2]) and we obtain 
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rribfv 



91 



mbfp : 
mbfs 



■mB[giEM - 



gh{E^' 
-MEb' 



Eb)]. 
-Eb)] , 



(4) 



where , -E'g' ^'^'^ Eb are the energies of the M meson, 
B' and B, respectively. In Eq. (|4]), 52 and /2 disappear, 
while (73 and /a are neglected since the corresponding 
terms are proportional to which is small comparing 
to itibEm and mB{E^, — Eb)- The form factors in Eq. 
([3]) can be related by the spin SU{2) and flavor SU{3) 
symmetries. In Table HI we show the relations for the 
form factors in {pp\{ub)v,A\B~''^). 



TABLE I: Relations between the B 
factors with i — 2, - ■ ■ ,5. 
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The decay width F of i? ppK is given by 



F = 



+ 1 p{niB-mK)' 
4mJ 



1/2 



inmB ) 



-\Afdtdcose, (5) 



where (3p = {1 - 4ml/ty/^, Xt 
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+ t' 



''K"^B' 



9 is the angle between the 
three-momenta of the K meson and the proton in the 
dibaryon rest frame, and is the squared ampli- 

tude of Eq. ([1]) by summing over all spins. Note that 
4tob£'p(p) = TUg + t — mp± PpX]^^ cos 9. From Eq. (O, 
we can study the partial decay width dT/d cos as a func- 
tion of cos9, i.e., the angular distribution. We may also 
define the angular asymmetry by 
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J— 1 dcosd 
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which is equal to {N+ — N-)/{N+ + N-), where N± are 
the events with cos 9 > and cos 6* < 0, respectively. The 
angular asymmetry in B~ ppK~ has been measured 
by the Belle Collaboration to be 



Ag{B ppK 



— U.OJ_Q 07 , 



(7) 



implying that the protons are emitted along the di- 
rection most of the time [2| in the pp rest frame, which 
seems to be unexpected from the previous B studies in 
the PQCD picture 

In general, when a decay mode mixing with vec- 
tor (axial-vector) and scalar (pseudoscalar) currents, it 
makes the partial decay width as a form of 



dV 
dcos 9 



~ a cos^ 9 + b cos 9 - 



(8) 



For B^ Ap7r+, one gets that c > — a > with 6 ~ 
[1^ . Therefore, its angular distribution as a function of 
cos9 is figured as a parabolic curve opening downward 
[21I . [2^ , which is consistent with the data ]2] . On the 
other hand, the angular distribution of B~ ppK~ Q 
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is gradually bent up as cos6' = — 1 shifts to +1 (see Fig 
3a in Ref. 0]), leading to an asymmetric Ag in Eq. ([7]), 
which is unexpected since the decay is dominated hyV-V 
and A ■ A contributions [l^, Clearly, the data for 

B~ —f ppK~ indicates a non-neglected cos 6 term which 
is comparable with the cos^ 9 one in Eq. ([5]), i.e., a ~ b > 
0. To find out a large cos term, it is important to note 
that the energy difference of the proton pair Ep — Ep is 
proportional to cos 6 and related to 175 and /s as seen from 
Eq. Q , which could provide new source of the angular 
dependence. We now summarize all possible cos 6 and 
cos^ 9 terms in ppK~ as follows: 

Vi ■ V5 , V4 ■ V5 , Ai ■ A4 , A4 ■ Ar, cx cos 9 , 
Vu5) ■ Vu5) , Au5) ■ A(5) , ^1 ■ ^5 , « cos^9, (9) 



where we have denoted the terms corresponding to gi 
and fi in Eq. ([2]) as Vi and Ai , respectively. The squared 
amplitude in Eq. (O is reduced to be 
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with (i, A^^) = (t, Xy^)/m%, {rhp,EK) = {mp, EK)/mB 
and 



K 



mb{ms + ruu) 



.(11) 



Here, to simplify our formula wc have ignored the current 
contribution of C{B^ ppK^) in Eq. ([1]), estimated 
around 1% of the measured one by u sing the values of 
the time-like baryonic form factors in [l7|. However, all 
terms are kept in our numerical calculation. 

In our numerical analysis, we take the 10 data points 
of dBr/dcos9 in B~ — *■ ppK~ from Fig. 3a in Ref. 2| 
(see also Fig. [Ij and measured decay branching ratios 2| 
of (5. 74 ±0.61) and (1.20 ±0.35) x lO"*^ for ppR- 
and 5° — > ppKs, respectively, as input values, along 
with the CKM matrix elements referred in Ref. [iot . 
(01,04,06) = (1.05,-0.0441-0.0072i,-0.0609-0.0072i) 
[23, m„(mb) = 3.2 MeV, ms{mb) = 90 MeV ^ and 
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FIG. 1: Branching fraction vs. cos 6 in the pp rest frame 
for (a) — > ppK^ (n^) with solid (dash) curve and (b) 
— > ppKs, where the data points (stars) in (a) and (b) are 
from Ref. 0. 



mb{mb) = 4.19 ± 0.05 GeV [SO!]. We note that the ratio 
of {(^k/oikY is around 15% which suppresses fv,s (or 
/i,4,5) terms in Eqs. ^ and ([3]). Although there is no 
surprise that Vi ■ Vi and Ai ■ Ai create cos^ 9 with an 
expected minus sign, the dominant contribution arises 
from the g| term with t around 4 — 25 GeV^ and as a 
consequence, pg and pg2 can be both positive with the 
same size and the condition of a ~ 6 > is fulfilled. 
By performing the fitting, iV|| p and M||'^ in Table U are 



determined to be 



iV|| = 127.1 ± 26.6 GeF^ N-^ 



-200.9 ± 51.9 G'eI^^ 



Af|| = -25.0 ± 15.4 GeF^, il/,^ = 227.3 ± 22.0 GeF'', (12) 

As seen in Fig. [1^. our result (solid curve) for 
dBr{B~ — > ppK^)/dcos9 as a function of cos 9 explains 
the data Q well and it is in agreement with the Belle 
data of Ag in Eq. ([7]). For B° ppKs, the angular 
distribution is shown in Fig. [T}j, which is also consistent 
with the data except the one point close to cos9 = 0.8. 
Clearly, more data at the ongoing B factories are needed. 
In addition, we predict 



AgiB- 
AgiB'' 



> ppiT ) 

ppKs) 



0.45 ±0.11 , 
-0.35 ±0.10 



(13) 



The result in Eq. for Ag{B^ ppir^) is antici- 

pated since the decay is similar to B^ ppK^ although 
it is dominated by the tree-diagram. It is interesting to 
point out that the minus sign for Ag[B'^ ppKs) in Eq. 
p3)) is different from the expectation of the fragmenta- 
tion mechanism [lH. We note that from our fitted form 
factors we obtain 



Br{B^ ppTT^) ^ (3.0 ± 0.4) X 10" 



(14) 



which supports the Belle measurement of (3.06^" '''^ 
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± 



0.37) X 10"^ [H, but is higher than the Babar data of 
(1.24±0.32±0.10)xl0-6 0|. We look forward to having 
the future experiments at Belle and BaBar to check our 
predictions in Eqs. ^3]) and p4|) . 

It is known that the asymmetric Dalitz plot for B^ 
ppK~ reported by the BaBar Collaboration 0, sup- 
ports the fragmentation mechanism [ll| but disfavors the 



4 





the National Science Council of Republic of China under 
the contract #s NSC-94-2112-M-007-(004,005). 



3 . 5 

(GbV) 



FIG. 2: Branching fraction of Br{B^ ppM~) / dm^M with 
{a) M = K and (b) M = n as functions of jtiba/ (B = p,p) 
with the dash, dot and solid curves representing (i) mpj^f- > 



(ii) 



< m. 



and (ii), respectively. 



pM- 



and (iii) difference between (i) 



gluonic resonance state [13, HI as well as other interme- 
diate state in the pole model [sij. We now examine if 
we can give a quantitative description on the Dalitz plot 
in the PQCD approach. To do this, we define the Dalitz 
plot asymmetry for B BB'Af as 



^ ^ r(»^B'M)> - r(TOBM)> 

r("^B'M)> + r(77iBM)> 



(15) 



where r(TOBM)> [r(?Tig/]y[)>] denotes the decay width 
for the range of rriBM > "^b'M ("^b'm > ™bm) divided 
by the line of »tibm = "^b'm ^^'^ Dalitz plot, ft is 
easy to see that for B~ — *■ ppK~, i.e., B = B', the 
Dalitz plot asymmetry is identical to the angular asym- 
metry in Eq. ^ due to the fact that these two asym- 
metries arise from the same source of pe in Eq. (|10[) 
with the relations of cos6' = /3j7^^~^^^(™^if- 



and t 



2ml 



A^i^U- ^'^PK-)- Explic- 
itly, in Fig. [2^, we show the decay branching fraction of 
ppK~ as functions of rripx- p^- with the dash, 
dot and solid curves representing (i) m^^- > rUp^- , (ii) 
TnpK- < '"T^pK- and (iii) difference between (i) and (ii), 
respectively. It is interesting to note that, as seen from 
the figure, the Dalitz plot asymmetry peaked around 4 
GeV is exactly the same as the data in Ref. Q . However, 
our prediction for B^ ^ ppK~ shown in Fi g. ^Bp is dif- 
ferent from the BaBar unpublished result in [20| like the 
decay branching fraction. Clearly, more precise data for 
the Dalitz plot distribution on the tt mode are needed. 

Finally, we remark that the form factors in Eq. ^ de- 
termined by the angular distribution in i?~ ^ ppK~ can 
be used to examine other experimental measured three- 
body baryonic B decays with a vector meson in the final 
state, such as i? ^ ppK* and B" KpJ/^ P,[3], which 
have not been theoretically explored yet. Furthermore, 
direct CP violation in B~ ppK^*^~ can be also inves- 
tigated. Moreover, our study on the angular distribution 
can be extend to the above modes as well as the decay 
of B~ Ap7, which has only been discussed in the pole 
model (3^ . 
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